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Abstract. Two dimensional condensed matter is realised in increasingly diverse 
forms that are accessible to experiment and of potential technological value. The 
properties of these systems are influenced by many length scales and reflect both 
generic physics and chemical detail. To unify their physical description is therefore 
a complex and important challenge. Here we investigate the distribution of 
experimentally estimated critical exponents, j3, that characterize the evolution 
of the order parameter through the ordering transition. The distribution is found 
to be bimodal and bounded within a window ~ 0.1 < /3 < 0.25, facts that 
are only in partial agreement with the established theory of critical phenomena. 
In particular, the bounded nature of the distribution is impossible to reconcile 
with existing theory for one of the major universality classes of two dimensional 
behaviour - the XY model with four fold crystal field - which predicts a spectrum 
of non-universal exponents bounded only from below. Through a combination 
of numerical and renormalization group arguments we resolve the contradiction 
between theory and experiment and demonstrate how the "universal window" 
for critical exponents observed in experiment arises from a competition between 
marginal operators. 
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1. Introduction 

New types of two-dimensional system on which meaningful physical experiments can 
be performed include optical lattices of trapped atomic gases pQ , magnetic surfaces [2] 
and "(5-doped" magnetic layers [3J. These add to a list of well established two 
dimensional systems that includes ultrathin magnetic films [4], atomic monolayers 
(both physi- and chemisorbed) [5j [6l [71 [8J , crystalline surfaces [9] , superconducting 
layers [TO] and arrays of interacting Josephson junctions [TT]. Recent theoretical 
developments on the concept of "extended universality" [T3], the effects of finite 
size P21 HI] , and the dipolar interaction [TBI HI] should be particularly relevant to 
understanding experiments on these systems, both old and new. 

The key experiment on two dimensional systems is to test the existence and 
temperature dependence of a magnetic or crystalline order parameter m(T). In cases 
where m can be measured experimentally (which excludes, for example, superfluid 
films [7J), this is invariably found to approximate a power law over a certain range of 
temperature: m ~ (T c — T)@ , where T c is the transition temperature. Theory predicts 
a limited number of possibilities for the value of the exponent (3, as dictated by the 
universality class of the system. In two dimensions crystal symmetries and consequent 
universality classes are relatively few. We show here that the Ising, XY and XY with 4- 
fold crystal field anisotropy (XY/14) are the three main experimentally relevant classes. 
The three and four state Potts models provide additional universality classes observed 
in experiments on adsorbed gaseous monolayers [TTJUS] and surface reconstruction [§]. 
For the Ising, three- and four-state Potts models, (3 — A, |, ^ respectively. For the 
XY model, one expects /3 — 0.23, a universal number that arises in the finite size 
scaling at the Kosterlitz-Thoulcss-Bcrezinskii (KTB) phase transition [19, 20 , though 
not a conventional critical exponent [21]. For XY/14, theory predicts a continuously 
variable critical exponent (3 cx l/ft.4 and thus a continuous spectrum of values when 
sampled over many real systems (see references [351 123] and this work) . 

We have tested these ideas by means of an extensive survey of experimental two 
dimensional critical exponents, including data for magnetic ultrathin films, layered 
magnets that exhibit a temperature regime of two dimensional behaviour [24l [25l [26l 
[27], order-disorder transitions in adsorbed gaseous monolayers [T7J [18], and surface 
reconstructions [9] . The results are presented in Figure Q] and in the appendix. As 
observed previously on more limited data sets [H dH [251 IIS] , the distribution of /3's is 
distinctly bimodal, with strong peaks at (3 — 0.12 and (3 = 0.23, as expected for the 
Ising and XY models. In several cases ideal Ising [24 ) [3D] and XY [25 ] I26 ] [27 ] [3T ] [32] 
behaviour has been confirmed in great detail by measuring thermodynamic quantities 
other than the magnetization. Likewise there is compelling evidence for Potts 
universality in several non magnetic systems [TTJ, |T51 [331 (Ml M ■ However the XY/14 
universality is more elusive. Of particular relevance to the present discussion are the 
ferromagnetic monolayer Fe/W(100) [32] and the layered ferromagnets RL^CrCLt [25] 
and K2CUF4 [26], easy plane systems which have been shown to exhibit the full range 
of ideal XY behaviour despite their 4-fold symmetry. Another very well characterised 
easy plane system with a 4- fold crystal field is the layered antiferromagnet K2FeF4 |31j . 
but this is not XY-like, with (3 — 0.15 intermediate between the XY and Ising 
values. Claims for XY/14 universality have been made for the ferromagnetic films 
Fe/[Au or Pd](100) [35], characterised only to a limited extent, as well as the order- 
disorder transitions of H/W(011) [36] and O/Mo(110) 37J, for which full sets of critical 
exponents are available. The behaviour of these candidates for XY/14 is seen to fall 
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Figure 1. Histogram of j3 values for all two-dimensional systems reported in 
Tables IA1I IA2I and IA3I The universal window is highlighted by the grey shading. 
Criteria for inclusion in the data set are discussed in the appendix. 



into two categories, which on closer inspection appears to be related to the strength of 
h^: those with weak h\ are XY-like with (3 « 0.23, while those with stronger hi have 
exponents in between the XY and Ising limits, 0.125 < (3 < 0.23. Most strikingly, there 
is no experimental evidence of the divergence of the exponent j3 implied by (3 cx I//14. 
Instead, most experimental data that cannot be ascribed to the Potts classes lies in 
a "universal window", bounded by the Ising and XY values. There are exceptions 
at the upper bound where crossover to three dimensional behaviour may increase the 
value of (3 upwards from 0.23 [211 EE]- However, it is clear from the histogram that 
the majority of systems are indeed encompassed in a limited range between the Ising 
and XY values. 



2. Universality Classes in Two Dimensions 

Before we address the main question of why the universal window exists, it is relevant 
to specify the occurrence and relationships between two dimensional universality 
classes. Considering first magnetic degrees of freedom, we ignore the possibility 
of truly Heisenbcrg behaviour, remarking that the broken translational symmetry 
inherent to layers or surfaces, combined with a condition of crystal periodicity, means 
all real systems have at least one p-fold axis, which necessarily introduces relevant 
perturbations. Thus, although pure Heisenberg behaviour may be observable over a 
restricted temperature range [39J, HQ] , it must give way to behaviour characteristic of 
the perturbations at temperatures near to the phase transition. These perturbations 
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Table 1. Classification of continuous transitions which can be observed in two- 
dimensional magnetic systems and in structural order-disorder transitions on 
surfaces. • indicates the occurrence of a particular universality class, whereas 
X indicates its absence. The special case of the square lattice dipolar system is 
discussed in the text. 



Universality Class 


Magnetic Systems Adsorbed Systems 


Ising 


• • 


XY 


• • 


XY/14 


• • 


3-state Potts 


x • 


4-state Potts 


x • 



take the form of axial anisotropy (either easy axis or easy plane) and p-fold in-plane 
anisotropy (p = 1,2,3,4 and 6). Easy-axis systems are generally Ising-like (despite the 
fact that the normal to the plane is usually a polar axis) while easy-plane systems 
with p — 2-6 should be described by the XYh p model. XY/12 is in the Ising class, 
whereas XY/i 3 is in the 3-state Potts class, although it is very unlikely in magnetic 
systems owing to time reversal symmetry (we found no examples). XY/14 constitutes 
a universality class distinct from the 4-state Potts class, while the phase transition in 
XYh§ is in the XY class [35]. Inclusion of the dipolar interaction on lattices other 
than the square lattice does not add extra universality classes. However the case of the 
square lattice must be regarded as an unsolved problem: perturbative calculations [41] 
and numerical results [42, 43, 44 suggest that the square lattice dipolar model belongs 
to XY/14, but the renormalisation group calculations of Maier and Schwabl indicate 
a different set of critical exponents [15] . The experimental data considered here are 
consistent with the former result rather than with the latter, but Maier and Schwabl's 
prediction could yet be born out on an as yet undiscovered ideal model dipolar system. 
At least as far as the existing experimental data set is concerned, we conclude that, 
for magnetic systems, there are only three main universality classes: Ising, XY and 
XY/14. 

The situation is essentially the same in non-magnetic systems [35J|3n] but with the 
additional possibility of the 3- or 4-state Potts classes due to competing interactions 
beyond nearest neighbour 46, 47J. Indeed, Schick [46] used arguments from Landau 
theory to classify the phase transitions of two dimensional adsorbed systems into only 
four classes: the Ising, XY/14, 3- and 4-state Potts. This set is supplemented by a 
chiral 3-state Potts class which shares conventional exponents with the pure 3-state 
Potts class [48] (hence for present purposes we shall treat these two cases as a single 
class). One result of the current work is that the pure XY class is also relevant to 
order-disorder transitions in adsorbed layers. Combining these observations we have 
five universality classes for structural systems and three for magnetic systems, as 
summarised in Table [T] 

3. Calculation of Critical Exponents 

The relationship between the Ising, XY, XYh p and clock models may be discussed 
with reference to the following Hamiltonian: 

Hp = - J X! cos (^i ~ fy) - h P X! C0S (P^)> (!) 

(iJ) i 
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in which the #j's are the orientations of classical spins of unit length situated on a 
square lattice with periodic boundary conditions and confined to the XY plane, J is 
the coupling constant and h p is the p-fold crystal field. It should be noted that unlike 
real systems, the lattice symmetry in computer simulations does not constrain the 
spin symmetry, and consequently the adoption of a square lattice does not restrict 
the generality of our arguments. In the limit h p — > oo, the Hamiltonian ([1]) is called 
a clock model, since Oi is restricted to discrete values evenly spaced around a circle: 
2%(n/p), n = 1, . . . ,p — 1. Jose, Kadanoff, Kirkpatrick and Nelson (JKKN) [22] 
have shown that for p > 4, h p is an irrelevant scaling field down to intermediate 
temperatures, with the result that fluctuations restore the continuous symmetry of 
the 2dXY model above a threshold temperature, leading to a KTB transition [20] 
and quasi-long range order over a finite range of temperature. Recently it has been 
shown [12] that a similar scenario remains valid even for infinitely strong crystal 
field strength, with the result that fluctuations restore continuous symmetry for p- 
state clock models with p > 4, although for 4 < p < 6 this occurs above the KTB 
temperature, Tkt- For p = 2 and 3, h p is relevant, leading to phase transitions in 
the Ising and 3-state Potts universality class respectively. I14, on the other hand, 
is a marginal perturbation [22 . A second order phase transition is predicted with 
non- universal critical exponents depending on the field strength. As /14 — > 00, XY/14 
crosses over to the 4-state clock model, which is equivalent to two perpendicular Ising 
models, and the transition falls into the Ising universality class [49] ■ The non- universal 
transition for XY/14 is hence bounded by the Ising universality class for large /14. 

The non-universal exponents of XY/14 can be calculated analytically within the 
framework proposed by JKKN. They showed that to describe the evolution of the 
KTB transition in the presence of a weak p-fold field it is sufficient to replace fl} by 
the generalized Villain Hamiltonian [22] [50] 



k B T ^ 



l--(0 i -0 j -2irm ij ] 



y ipn$i 



+ log(yo)^2S 2 R + log(y p )Y^nl (2) 

i i 

where K = J/fc^T. The integers m%j maintain the periodicity of the original 
Hamiltonian, for rotations over an angle 2tt. Sr is a directed sum of integers my 
around a square plaquette of four sites centred at R: Sr = 77141 + mx2 — 17132 — 77743) 
takes values, Sr — 0, ±1, ±2 . . . and is therefore a quantum number for a vortex of spin 
circulation centred on the dual lattice site R. yo is related to the chemical potential fi 
and fugacity y for the creation of a vortex anti- vortex pair on neighbouring dual lattice 
sites: y = yo exp(— flu) ss yo exp(— tt 2 K/2). In the original Villain model yo = 1 but it 
is introduced here as a phenomenological small parameter which is renormalized in the 
subsequent flows. Similarly y p is a fugacity for a locking process of spins along one of 
the p-fold field directions with integer being a measure of the strength of this process 
at site i. For weak crystal fields, y p = \h p with h p = (h p / kgT), which reproduces 
the field contribution to the partition function to leading order in y p . For strong 
fields y p — > 1 and ([2]) transforms into a discrete p-state model. Note however that 
this is not the p-state clock model: although the Villain model maintains the global 
rotational symmetry it destroys the local 0(2) symmetry of the pair interaction. The 
discrete terms (8i — 9j — 27rmy ) 2 , Bi — (n/p)2ir : n — 0, 1, ... ,p — 1 hence do not 
have this symmetry over the interval — ir < (8i — 8j — 2 / nmij) < tt. For p — 4 this 
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means that neighbouring spins orientated perpendicularly have an energy less than 
half that of antiparallel spins and the ordered state has lower lying excitations than 
the corresponding clock model. It is therefore not clear whether the Villain model 
falls into the correct universality class in the strong field limit and for quantitative 
studies one should use Hamiltonian flj rather than ([2]). 

With yp set equal to zero, a direct space renormalization analysis for the spin- 
spin correlation functions resulting from (|2|) leads to RG flow equations for an effective 
coupling constant K c s and vortex fugacity, y. For K e g = 2/tt, y — 0, the flows yield 
the KTB transition [ST]. In the presence of the p-fold field the flow equations are 
modified and a third equation is generated [22l[52]. For the explicit case with p = 4, 
these are 

(X" 1 )' = K~ l + 4 {7T 3 y 2 e-^ K - A^R^yle^") ln(6) (3a) 
y'o = y + (2-7rK)y ln(b) (36) 
2/4 =2/4+ (2-— — Jy 4 ln(6), (3c) 

where b is the scale factor and where the equations are valid as b — > 1. This 
set of equations has fixed points at K* = 2/ir, j/q = ±2/1- We can calculate the 
linearized transformation matrix evaluated at the fixed point, *: M^j = fj^l*, where 

Ki = K~ 1 ,yo,y4 : - 

Solving for the eigenvalues we find 

A = l,l + f ±^V4a 2 + a 2 , (4) 

where a = 16ir 2 (2ir - l)y 2 e" 27r ln(6), a 2 = 2-yd, 6 = fyln(&), 7 = 8ir 3 ye- 27r ln(6), and 
where j/o — J/4 — V- Writing A = b a we extract the three scaling exponents. There 
is one relevant exponent, which is interpreted as <r\ = 1/V, the exponent taking the 
coupling constant away from the critical value at the now regular second order phase 
transition. There is also one irrelevant variable 02, which is interpreted as driving 
the vortex fugacity to zero. Finally, there is one marginal variable, 173, which, as 
announced, corresponds to the scaling exponent of the 4-fold crystal field. Taking /14 = 
all eigenvalues become marginal, consistent with the particular scaling properties of 
the 2dXY model. In the small field limit, o\ = —0-3 = 47re~ 7r /i4 and 02 = 0. This gives 
the non-universal correlation length exponent [22] v w I.S^bTkt/^)- The strong 
field limit, 2/4 = 1 gives v w 0.47, which should be compared with the exact result 
for the Ising model, v = 1. The agreement is poor, as might be expected given the 
distortion of the four fold interaction imposed by the Villain model. It is clear from 
this result that a quantitative calculation for the strong field limit requires a different 
starting Hamiltonian. 

In order to calculate (3 from the scaling relations [53], a second relevant scaling 
exponent is required. In this case the anomalous dimension exponent 77 can be 
calculated directly from the correlation function [22] . At the KTB transition of the XY 
model, rj — 1/4, giving the universal jump in the effective spin stiffness, K g — 2/tt. It 
follows from the scaling relation 2(3 = (d — 2 + rj)v that the finite size scaling exponent 
0/v = 1/8, as in the Ising model, despite the fact that here the true (3 and v are not 
defined. This is an example of "weak universality" [S3] between the two models. A 
striking result in the presence of a 4-fold field is that 77 remains unchanged to lowest 
order in /14 [52] , indicating that a weak universal line extends out from the XY model 
along the /14 axis. Here we make the hypothesis that the line extends right to the Ising 
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limit, in which case 77 = 1/4 for all h 4 . This is clearly a reasonable assumption for the 
level of calculation made here. It is also an appealing result as other examples of weak 
universality are far less accessible to experiment [55] . Analysis of the numerical data 
presented in the next section lends weight to this hypothesis, although the observed 
behaviour is found to divide into two regimes, depending on the strength of the /14 
field. 

From this analysis we therefore predict a range of non-universal magnetization 
exponents going from 



8 V hi 

for weak field, to [3 — 1/8 in the strong field limit. To make quantitative comparison 
with simulation and experiment we need to estimate (3 as a function of /14/J. The 
critical value K* — 2/ it corresponds to a renormalized coupling constant, J c ff, valid 
at large length scale such that fc^TxT = In general J c s < J: for the Villain 

model fc^TKT ~ 1.35J [56] , while for the XY Hamiltonian |lj fesTxT/ 'J ~ 0.9 and is 
different again for more realistic Hamiltonians. Hence, while we can make a theoretical 
prediction for the low field behaviour, 



/14 

with a a constant of order unity, scaling equation j5| by a factor JcbTkt/J will 
probably not lead to an accurate quantitative estimate for a and the precise value 
is beyond the scope of the present calculation. 

4. Competition with Essential Finite Size Effects 

The survey of the (3 values illustrated in Figure [1] shows a clear discrepancy between 
theory presented above and experiment: the large values of (3 predicted for small /14 
do not appear and the range of values is cut off at f3 « 0.23. As the latter is an 
effective exponent characteristic of XY criticality up to a finite length scale, it seems 
clear that the non-universal critical phenomena are suppressed, for weak field, by the 
exceptional finite size scaling properties of the pure 2dXY model [3TJ [571 [55] . This 
hypothesis can be tested by numerical simulation, in which both /14 and the system 
size may be directly controlled. 

In a real XY system the relevant length scale will in most cases be less than the 
physical size of the system: for example, it could be a coherence length controlled 
by defects or dipolar interactions [59j . or, in the case of layered systems, a crossover 
scale to the third dimension [501 [2T]. Thus, although real systems might have, for 
example, 10 16 spins, the relevant scale for XY critical behaviour will typically be 
much smaller and compatible with the scale of Monte Carlo simulations, where the 
appropriate length scale is simply the system size. This finite length scale gives 
rise to a finite magnetization that disappears at the rounded KTB transition. As 
emphasised in [3T] , this is perfectly consistent with the Mermin- Wagner theorem [FT] , 
which proves that the magnetization will be strictly zero in the thermodynamic limit. 
It is easy to convince oneself that finite size corrections to the thermodynamic limit 
are important for any physically realizable cut off length scale. The resulting low 
temperature magnetization is therefore directly relevant for experiment. 



In Figure 2(a) we show the magnetic order parameter, the thermally averaged 



magnetic moment normalized to unity, versus temperature, with different 4-fold field 
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Figure 2. Monte Carlo data for the 2dXY model in the presence of a 4-fold 
crystal field, with N = 10 4 . Plot (a) displays the magnetization data of the XY 
model in the presence of increasingly strong anisotropics. Plot (b) displays the 
same data as a function of reduced temperature t, in a logarithmic scale. In both 
plots the "universal window" is highlighted by the blue shading. 



perturbations, for a system with N = 10 spins. For h± — the magnetization is 
characterized by the effective critical exponent, (3 ph 0.23. A finite size analysis of 
Kosterlitz' renormalization group equations shows that in the region of the transition 
it approaches a universal number (3 = 37r 2 /128 ~ 0.23, in agreement with both 
experiment and simulation data, such as that shown here. For weak crystal field, 
hi there is no change in the the region of the transition and the magnetization 
data coincide with the data for zero field [62]. Only for J14/J > 0.5 do they leave 
the zero field data through the transition, approaching results for the 4-state clock 
model for large values of /14/J. In Figure [2(b)] we show log(m) against log(i), where 
t = (T — T c ) /T c . The transition temperature T c is calculated from a finite size scaling 
analysis of the fourth order Binder cumulant for M [53] [MJ [55] and is an estimate of 
the value in the thermodynamic limit. The slopes, for small t, give a first estimate of 
the exponent (3, indicating that it lies in the interval 1/8 < j3(h±) < 0.23 for all values 
of /14, exactly as observed in experiment. The crossover to Ising behaviour is slow: for 
/14/J = 1, /3Q14) w 0.15 and to approach (3 w 1/8 requires /14/J in excess of 5. 




Figure 3. Best data collapses for the 2dXY model with 4-fold crystal field 
h^/J = 1.0 and h^/J = 2.0 for different system sizes. 



Hence the data here, as in previous numerical work [SH [551 IM] , show evidence for 
a finite pocket of XY critical behaviour for small values of hi. This appears to refute 
the prediction of JKKN, derived explicitly in the previous section, that the exponents 
vary continuously with hi [65] (see the further discussion below). For intermediate 
field strengths, however, the non-universal criticality does appear to hold as can be 
confirmed by a more detailed finite size scaling analysis. The values of j3 and v can be 
estimated more accurately by collapsing data for various system sizes onto the scaling 
relation ML"/" ~ tL 1 ^. The best data collapses for hi/ J = 1 and hi/ J = 2, with T c 
in each case fixed from the Binder cumulant calculation, are shown in Figure [3J We 
find j3 — 0.148(5) and (3 = 0.136(10), in good agreement with the values found from 
Figure |2(b)[ and v = 1.19(4), v = 1.09(8). The ratio (3/v = 0.126(4) in each case, is in 
agreement with the weak universality hypothesis. Similar results for hi/ J — 0.5 can 
be found in [65) . Although these exponent values are not so different from the Ising 
model values, the data collapse is less satisfactory when Ising exponents are used. 

Further evidence for weak universality at intermediate field strengths can be found 
from studying the finite size scaling properties of m at the transition. In Figure [4] 
we show log(m) against log(L) for hi/ J = 1 for a range of temperatures near the 
transition. At the transition one expects a power law evolution characterised by the 
finite size scaling exponent r)/2 — j3/v. The best power law occurs at T c — 1.010(5) J, 
which is the same as the value found from the Binder cumulant method. The scaling 
exponent rj/2 = 0.126(3), is the same as that found for the data collapse in Figures 
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Figure 4. Magnetisation m against system size L, in log-log scale, for the 2dXY 
model with 4-fold crystal field /14/J = 1.0, for temperatures near the transition. 



3(a)] and [3(b)] 



This and previous numerical work [5H [551 [531 [55J are consistent with /14 being 
marginal. In this case the crossover exponent to the new universality class is zero so 
crossover occurs, at best on exponentially large length scales, as a result of corrections 
to scaling [62] • Hence, for small and intermediate crystal field strengths the finite 
size scaling appears compatible with that of the continuous symmetry of the 2dXY 
model [S^, as in the 6-fold case. In fact the most detailed finite size scaling analysis [53] 
shows no evidence of such a crossover for small fields. It therefore remains an open 
question whether the pocket of pure XY behaviour for small h±/J is a pragmatic 
observation related to excessively slow crossover, or whether it remains right to the 
thermodynamic limit. In either case this is the main result of this section: large values 
of f3 are indeed masked by the pocket of 2dXY behaviour, leading to the effective 
exponent (3 for weak /14 and creating a divide between systems with strong and weak 
4-fold fields, with the non-universal character of XY/14 only appearing for f3(hi) < 0.23. 
The threshold value of /14, separating the two regimes can be estimated theoretically 
by putting — 0.23 in Equation ([6]). Using ksT^/J sa 0.9 gives a = 1.6 and 

hi/ J ~ 0.9, a ratio of order unity, in agreement with the above general arguments, 
but an over estimate compared with numerics, where the change of regime occurs for 
hi/ J ~ 0.5, corresponding toawl. 

Having confirmed that 77 w 0.25 over the whole range of /14, we finally fix rj = 0.25 
and use our estimates of v(h±) from the scaling collapse to give a further estimate of the 
exponents as a function of /14. The estimates of v and /3, summarised in Tabled are in 
good agreement with all previous unconstrained estimates. We also include estimates 
of derived by a typical experimental analysis of fixing T c from the maximum in 
the susceptibility or where the magnetization approaches zero, and deriving (5 from 
a log-log plot. There is seen to be a systematic error between the different estimates 
of P, especially for small values of /14. Nevertheless, the experimental exponents are 
still found to lie in the universal window of values predicted for the "true" exponents 
of the underlying model. The various critical exponents plotted in Figure [5] are found 
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Table 2. Critical exponents for the XY/14 model, as determined from a finite size 
scaling analysis, and as measured directly from Monte Carlo magnetization data 
for a system of size L = 100. 

hi v (3 T C L f3(T c L ) 

0.5 1.37(6) 0.171(10) 1.01(1) 0.214(9) 

1 1.19(4) 0.148(5) 1.04(1) 0.196(6) 

2 1.09(8) 0.136(10) 1.08(1) 0.155(3) 
5 1.04(6) 0.130(7) 1-12(1) 0.129(3) 
00 1.00(5) 0.125(6) 1-14(1) 0.123(3) 




to be linear in I//14. By fitting to j3{hi) = 0.125 + a/h^, we estimate the constant a 
to be 0.032 for the "true" exponents, and 0.05 for the experimental exponents. These 
values are clearly very different from that expected for the constant a in Equation ([6]) , 
but once outside the pocket of pure XY behaviour we are no longer in the weak field 
regime for which Equation ^ is valid, as was shown in the previous section. 

5. The Strong Field-Weak Field Divide 

Experimental evidence for the strong field-weak field divide comes from making head 
to head comparisons between systems listed in the appendix. 

A quantitative comparison is afforded by the ferromagnet Ri^CrCLj [57] and 
the antiferromagnet K2FeF 4 [3T], both quasi two dimensional square lattice systems 
with S = 2. In both systems the intra-plane isotropic exchange coupling, J is much 
bigger than the inter-plane value J', giving rise to an extended temperature range 
with two dimensional critical fluctuations. However, while the ferromagnet shows all 
the characteristics of the pure XY universality class [35J HB] , the antiferromagnet has 
non-universal exponents, with /3 = 0.15 [31] , which we now see to be consistent with 
XY/14. A realistic model Hamiltonian for either system has the following form 

n = Jo s * ■ s i + ^E(^) 2 + \< S X + 5 -)> ( ? ) 
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Table 3. Main parameters for K2FeF4 and Rb2CrCl4, as determined from 
experiment 1311 1671 . 





K 2 FeF 4 


Rb 2 CrCl 4 


s 


2 


2 


J(K) 


-15.7 


15.12 


D(K) 


5.7 


1.06 


E(K) 


-0.49 


0.123 


D 


0.363 


0.07 


e 


0.0052 


0.0013 


P 


0.15(1) 


0.230(2) 



where the weak interplane exchange and a weak fourth order axial term are ignored 
(in the case of Rt^CrCU small departures from tetragonal symmetry are neglected for 
the purpose of this discussion) . 

The crystal field D confines the spins to an easy plane breaking the 0(3) rotational 
symmetry of the Heisenberg exchange and the 4-fold term e breaks symmetry within 
that plane, putting Hamiltonian |(SJ) in the XYft, 4 universality class. For both lU^CrCU 
and K2FeF 4 accurate estimates of the Hamiltonian parameters were derived by fitting 
magnon dispersions measured by neutron scattering to a self-consistent spin wave 
calculation [3TJ [BTl • However, in order to fit the spectra, the fourth order term was 
decoupled into an effective second order term, with amplitude E m 6eS 2 . In the low 
temperature limit one can estimate parameters e = |eS' 4 /Jo5' 2 | and D = \DS 2 / JqS 2 | . 
Values are shown in Table [3] for both materials. 

To get an estimate of the 4-fold field that determines the critical exponents, 
it is tempting to assume that systems with S — 2 are classical and to associate e 
with the parameter hi/ J arising from Equation fl}. The parameter e = 0.0013 for 
Rb2CrCi4 and 0.0052 for K2FeF4, which seem sufficiently small to put both systems 
into the weak field regime with pure XY universality, as has been directly confirmed by 
numerical simulation [70] . However, assigning an effective classical Hamiltonian of the 
form ([1]) for systems with finite S is not so straightforward: for finite S, through the 
uncertainty principle, out-of-plane and in-plane spin fluctuations are not statistically 
independent. As a consequence the energy scale for in-plane spin rotations and the 
consequent effective value for hi depend collectively on J and D as well as on e. 
This can be seen from a detailed consideration of the magnon dispersion arising 
from iJSJ. This calculation reveals a distinct difference between the ferromagnetic 
and antiferromagnetic cases, with the latter retaining strong quantum effects even for 
S = 2. For the antiferromagnet one finds two magnon branches which, for D = e = 0, 
are degenerate and gapless for zero wavevector and where each mode constitutes a 
conjugate in-plane and out-of-plane spin fluctuation term of equal amplitude. For 
finite crystal field strength the degeneracy is lifted, energy gaps appear everywhere 
in the spectrum and the symmetry is broken between the in-plane and out-of-plane 
fluctuation amplitudes. In the following we refer to a mode as in-plane or out-of-plane 
if the conjugate variable with the largest amplitude is in, or out of the plane. To 
lowest order in 1/S the out-of-plane branch develops a gap at zero wave vector: 

A x = S [2(D + \E\)(2\J\z + 4\E\)} 1/2 » 2Sy/Dz\J\, (8) 

while the in-plane branch has 

A 2 = S [(4|£|)(2| J\z + 2D + 2\E\)} 1/2 « 2Sy/2\E\z\J\, (9) 
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(here z — A). These gaps depend on the geometric mean of the exchange field zJ 
and the crystal field e or D with the result that they are surprisingly large on the 
scale of J, as noticed by Thurlings et al. [31] • For K 2 FeF 4 Ax = 70.8 K, A 2 = 23.9 
K at 4.2 K, renormalising only weakly with temperature |31j . Ai is larger than 
the transition temperature, Tn = 63 K ss JS 2 , so the out-of- plane branch of spin 
fluctuations will be frozen by quantum effects over the whole of the ordered phase, 
leaving the predominantly in-plane spin fluctuations only. Interpreting these as the 
classical fluctuations in an effective plane rotator model with Hamiltonian fl} leads 
to a crystal field, h^eS), of the order of A 2 . This gives fi4(e&)/JS 2 ~ 0.33 which 
is the right order of magnitude to fall into the strong field category. The fact that 
for K 2 FeF4, T^/JS 2 « 1, as is the case for the model systems with Hamiltonian 
(QJ presented in the previous section, is highly consistent with this interpretation. 
Higher order terms in l/S renormalize D and \E\ such that the values given in table 
[3] are higher than those predicted by fitting to linear spin wave theory. One can 
further speculate that quantum fluctuations for the in-plane branch will renormalize 
the effective /14 in (JT|) [64] to an even higher value. The non-universal exponents 
observed for K 2 FeF4 could therefore be examples of the XY/14 universality class. 

For the ferromagnet Rb 2 CrCl4, D flattens the cone of spin precession giving 
a range of q values where the energy spectrum varies approximately linearly with 
wave vector, but does not open a gap. The field e opens a zero wave vector energy 
gap that varies as V De. It is of order 1 K, decreasing to zero at the transition 
temperature T c — 52 K, and so can hardly affect the thermodynamics in the critical 
region. Although the effective value of /14 depends on the geometric mean of D and e 
rather than just the bare value of e it is independent of J and hence much smaller than 
than for the antiferromagnetic case. This places Rb 2 CrCl4 in the weak field regime, 
consistent with the observation of XY universality for this material [25] [68] . 

From this comparison, it seems likely that magnetic systems that show true XY/14 
universality will mostly be antiferromagnetic. Indeed a similar "spin dimensional 
reduction" due to quantum suppression of fluctuations has recently been observed 
in quantum Monte Carlo simulations with a Hamiltonian similar to ^ [58) . More 
calculations beyond the spin wave approximation are required to clarify this point. 

Among non-magnetic systems, oxygen absorbed onto Mo(110) [37] or W(110) [71] 
and hydrogen on W(110) 36] have both been claimed to fall into the XY/14 class, 
representing four fold equivalents of the two stage melting process for hexagonally 
coordinated systems [45]. Note that the (110) surface does not have four fold 
symmetry, but if we adopt these claims as a premise, then a comparison of the two 
systems is indeed perfectly consistent with XY /14 universality and with the preceding 
arguments about the strong field-weak field divide. Electron hybridization between 
absorbed and substrate particles will result in the generation of electronic dipoles 
aligned perpendicularly to the (110) surface. The resulting 1/r 3 interaction between 
the particles is repulsive and of sufficiently long range to ensure crystalization into a 
square lattice. The (110) surface provides a substrate potential with 4-fold topology 
(though not 4-fold symmetry) and which can be made commensurate with the free 
standing array by tuning the adsorbate density, the clearest example being the (2 x 2) 
lattice structure [36]. The result is claimed to be in the XY/14 universality class and 
the measured exponents, (3 ~ 0.19 [37] [71], are, in light of the current work, consistent 
with this. In principle, the same should be true for the (2 x 2) ordering transition 
for hydrogen on W(110) but the measured /?, 0.25, is consistent with the pure XY 
model [36] . Hydrogen being so much lighter that oxygen, larger zero point fluctuations 



Universal Window for Two Dimensional Critical Exponents 



14 



should make the substrate potential less effective at pinning the crystal, putting it 
in the category of systems with a weak field /14, consistent with the experimental 
observation. 

6. Other Exponents and Scaling Relations 

Further evidence for the experimental relevance of the finite size effects is found in 
the behaviour of other critical exponents. The exponent 77, which according to the 
previous weak universality arguments should be 0.25 for all /14, is only found to closely 
approximate the theoretical value for model Ising systems such as Rb2CoF4 [72]. 
For model XY systems the predicted 77 — 0.25 or S = 15 are always observed 
at a temperature well below T c (say 0.9 T c ), with rj{T) increasing to larger values 
at T c and S(T) decreasing, since S = (4 — r\)jr\. For example in the XY layered 
ferromagnets Itf^CrCLi and K2CUF4 77 (T) and 5(T) have been measured with precision 
by several different methods [23 [731 EH] : m both cases 77 rises to about 77 = 0.35 
at T c . This is a very strong signature of the finite size scaling properties of the XY 
model and is consistent with the predicted logarithmic shift in transition temperature, 
[T C (L) — Tkt] ~ iog^(L) ED EI]' f° r a nn ite size system [74]. As the measured value 
of 77 increases continuously through the transition its value at T C (L) is thus expected 
to be in excess of 77 = 1/4. 

It seems that the anomalous value of rj > 0.25 extends to systems with XY/14 
universality: for example, in K2FeF4 it is estimated to be 77 ss 0.35. This is again 
consistent with the shift in transition temperature observed in finite size systems. 
Defining T C (L) from the maximum susceptibility leads to a shift, [T C (L) — T c ] ~ L~ l ' v . 
Here, in the four fold field problem, v > 1 which means that shift remains important 
even in the intermediate field regime. Referring to Figure^ one can see that extracting 
a critical exponent from the initial slope, for T > T c , will lead to an overestimate of 
77. As experiments do not, in general, have access to the finite size scaling information 
available to numerical studies, it seems reasonable that the experimental 77 values are 
generally larger than the expected thermodynamic limit value. Thus, we propose that 
Tj(T c ) appearing greater that 1/4 remains a finite size effect. 

Similarly, the measured values of v are systematically smaller than unity, while 
a consequence of weak universality is that v should be greater than one for all finite 
/14. For example for K2FeF4 v as 0.9, giving j3/v = 0.16, greater than the predicted 
ratio 1/8, but together with 7 as 1.5 the set of exponents do satisfy the hyperscaling 
relation, 2/3 + 7 = dv, as well as the relation f3/v = d — 2 + 77/2. The same holds true 
for oxygen on W(110) [7T], for which f3 and 7 have been determined to be 0.19 and 
1.48, respectively. The shift in v is therefore consistent with the shift in 77. It seems 
reasonable to assume that these changes are also due to finite size effects, which at 
present prevent the observation of the weak universal line we have shown evidence for, 
for all values of /14. More detailed experimental and numerical studies to clarify this 
point would be of great interest. 

7. Conclusions 

In conclusion, the XY model with four fold crystal field is of relevance to a 
great number of experimental two dimensional systems. We have focused on the 
largest experimental data sets, those for two-dimensional magnets, adsorbed gaseous 
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monolayers and in particular on the measured exponent (3. With regard to the 
histogram in Figure [TJ the systems that comprise it can only be fully understood 
on a case by case basis. 

However, we show in this paper that the Hamiltonian (TTJ) contains the principle 
two dimensional universality classes that are relevant to experiment and that a uniform 
distribution of values /14 would, because of the marginal finite size scaling properties 
of the model, produce a probability density of the same form as Figure [1] with a 
continuous spectrum bounded by peaks at the Ising and XY limits. This is what we 
refer to as the universal window for critical exponents. We have further shown that the 
actual values of the four fold crystal field that occur in real systems are, at first sight, 
too small to take any system away from the XY limit. However, we have identified 
at least one mechanism, in antiferromagnets, whereby the four fold field is effectively 
amplified by quantum confinement of the spins to the easy plane. Other mechanisms 
of realising XY/14 universality are possible in individual cases [3TJ [37]. We have 
demonstrated the relevance of finite size scaling corrections to the experimental data 
set, with the relevant length scale giving a crossover away from XY criticality. Future 
work should focus on the finite size scaling aspects and on individual systems to see 
if a more accurate quantitative connection between the physical /14 and the observed 
critical behaviour can be established. Further to this, we propose here that the non- 
universal exponents of XY/14 should satisfy weak universality, with (3/v = 1/8 for all 
hi and we have given evidence that this is true in the range of intermediate field values. 
The robustness of the pocket of true XY behaviour, observed for weak fields [62|,l65]. 
in the thermodynamic limit remains an open question. Finally we remark that all 
evidence confirms that truly two dimensional systems, quasi two dimensional systems 
and numerical simulations reveal the same syndrome of behaviour, so much can be 
learned about new two dimensional systems [2l [3] through comparisons with old 
ones }311 167] . It is fortunate that there is such an extensive and carefully determined 
data base. 
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Appendix: Construction of the Histogram of [3 Exponents 

In constructing the histogram of experimental two dimensional [3 exponents, a number 
of factors were considered. First, it was crucial to avoid circular logic by excluding 
those systems which were assigned a dimensionality purely on the basis of their 
exponent values, rather than on a large body of experimental evidence. Fortunately, we 
found no such cases in the literature. Therefore all systems included in the histogram 
are assigned as two dimensional on the basis of compelling experimental evidence of 
two dimensionality Likewise we found no examples of systems considered to be two 
dimensional that exhibit (3 « 1/3, which might, in the absence of extra evidence, 
be mistakenly assigned as three dimensional systems and wrongly excluded from 
the data set. It should be noted that in layered magnets, the crossover from two 
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dimensional to three dimensional exponents is generally very sharply defined so there 
is no ambiguity in identifying the two dimensional regime. A second criterion for 
inclusion in the histogram was that the exponents were determined with reasonable 
precision and accuracy (typically A/3 < ±0.01). This inevitably necessitated a 
subjective judgement, but only a few results were excluded on these grounds: typically 
those exponents determined by powder (rather than single crystal) neutron diffraction, 
which is generally accepted to be inadequate for the accurate determination of 0. The 
experimental exponents are generally not asymptotic exponents, but the numerical 
study presented above reveals that the difference between asymptotic exponents 
and those determined using finite size scaling techniques at temperatures down to 
~ 0.9 T C L is generally negligible at the level of accuracy required for the present 
purpose. The histogram also excludes a number of interesting systems on the basis 
of there being legitimate grounds for alternative explanations for their observed 
critical behaviour. These include metamagnetic materials [75], systems undergoing 
spin-Peierls transitions [75] [77J [75], and bulk systems undergoing order-disorder 
transitions [79] . 

The following tables lists all the systems included in the histogram. Table IA1I 
contains data for layered magnets, and includes examples of molecular magnets (80;. 
Note that K2MnF4 represents two data points in the histogram as the elegant work of 
van de Kamp et al. [HI] used a magnetic field to tune the system between Ising and XY 
symmetry, with (3 recorded for both cases. In all other cases the /3's are determined 
in zero applied field. Table IA2I contains data for ultrathin magnetic films. Although 
there are several cases in which films of different thicknesses have been measured in 
order to study crossover to three dimensional behaviour, only the values of (3 in the 
two dimensional limit are reported here, and are included as only one data point in 
the histogram. Finally, Table lA"3l includes data for adsorbed gaseous monolayers, and 
systems undergoing surface reconstruction and melting processes. 

Our aim has been to make an exhaustive survey up to the time of publication. We 
apologise to any authors whose work we may have inadvertently overlooked, but we are 
confident that these cases would not significantly modify the form of the histogram. 
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Table Al. List of 2d critical exponents (3 for layered magnets reported in 
the literature, mostly measured by neutron diffraction (F = ferromagnet, A 
= antiferromagnet, Fo = (HCO2), chdc = trans-l,4-cyclohexanedicarboxylate, 
5CAP = 2-amino-5-chloropyridinium). 



System 


P 




i range 




T c (K) 


Type 


Reference 


Rb 2 CoF 4 


0.119(8) 


1 


1 n—l 

10 — 


2- 


10" 


4 


102.96 


A 


1821 


ErBa 2 Cu307 


0.122(4) 




nil 1 
J. 11 — 1 


■ 10- 4 




0.618 


A 


83 

. . 


K2C0F4 


0.123(8) 


1 


1 n — l 
10 — 


8 ■ 


10 


-4 


107.85 


A 


24 


BaNi 2 (As0 4 ) 2 


0.135 


3 


1 n-i 
1U — 


1 




-2 


19.2 


A 


27 


Ba 2 FeF 6 * 


0.135(3) 


7 


1 n—l 
10 — 


A 

1 ' 


in- 
1U 


-3 


47.96(4) 


A 


84" 


K 2 NiF 4 


0.138(4) 


2 


1 n-i 
1U — 


Q . 

O 




4 


97.23 


A 


85 


K3Mn 2 F7 


0.154(6) 


1 


1 n-i 
1U — 


1 

1 ■ 


IU 


-3 


58.3(2) 


A 


[86] 


Rb 2 MnCl 4 (B < 5.8 T) 


0.15(1) 


1 


-LU 


J . 




-3 


54 


A 


EE] ig 


K 2 MnF 4 


0.15(1) 


1 


1 n-i 

1U — 


1 . 


10" 


-3 


42.14 


A 


[88] 


K 2 FeF 4 


0.15(1) 












63.0(3) 


A 


31 


Rb 2 MnF 4 


0.16(2) 


1 


10" 1 - 


3 . 


10" 


3 


38.4 


A 


85 


Pb 2 Sr 2 TbCu 3 8 


0.165(5) 












5.30(2) 


A 


89 


BaFeF 4 


0.17 


3 


10" 1 - 


1 . 


10" 


2 




A 


90 


Cr 2 Si 2 Te 6 


0.17(1) 


6 


lO" 1 - 


3 . 


10" 


-2 


32.1(1) 


F 


91 


CsDy(Mo0 4 ) 2 


0.17(1) 












1.36 


A 


92 


0oCl 2 -6H 2 O 


0.18 


4 


lO" 1 - 


4 ■ 


10" 


2 


2.29 


A 


90 


MnC 3 H 7 P03-H 2 0t 


0.18(1) 


4 


10" 1 - 


1 . 


10" 


-2 


~ 15 


F 


93 


MnC 4 H 9 P0 3 ■H 2 O t 


0.18(1) 


4 


-LU 


2 • 


10" 


2 


~ 15 


F 


[93] 


KFeF 4 


0.185(5) 


3 


in- 1 

-LU 


j . 


10" 


-2 


137.2(1) 


A 


[92 [90] 


Fe(NCS) 2 (pyrazine) 2 


0.19(2) 


2 


10- 1 - 


3 . 


10" 


2 


6.8 


A 


[95] 


Rb 2 FeF 4 


0.2 


3 


10-1- 


2 • 


10" 


3 


56.3 


A 


85 


La 2 Co0 4 


0.20(2) 












274.7(6) 


A 


96 


MnC 2 H 5 P0 3 -H 2 0t 


0.21(2) 


6 


lO" 1 - 


9 . 


10" 


-2 


~ 15 


A 


[93] 


NH 4 MnP0 4 ■ H 2 Ot 


0.21(3) 


8 


10-1- 


2 • 


10" 


-2 


17.5(1) 


A 


[93] [97] 


K 2 CuF 4 


0.22 


3 


10-1- 


3 . 


10" 


-2 


6.25 


F 


[98] 


CuFo 2 -4D 2 Ot 


0.22(2) 


5 


10" 1 - 


5 . 


10" 


-2 


16.72 


A 


99 


CuFo 2 -CO(ND 2 ) 2 -2D 2 Ot 


0.22(1) 


4 


10-1- 


1 . 


10" 


-3 


15.31 


A 


[99] 


Tanol suberate§ 


0.22 


7 


1Q-1- 


2 . 


io~ 


-2 


0.7 


F 


100 


Sr 2 Cu0 2 Cl 2 


0.22(1) 


2 


io-i 

-LU 


1 . 


10" 


-2 


265.5(5) 


A 


[40] 


MnFo 2 • 2H 2 


0.22(1) 


4 


io-i 

-LU 


4 . 


io~ 


-2 


3.6 


A 


101 


La 2 Ni0 4 


0.22(2) 


8 


10 — 2 - 


2 . 




3 


327.5(5) 


A 


102 


BaNi 2 (P0 4 ) 2 


0.23 


3 


10 _1 - 


2 . 




2 


23.5(5) 


A 


[27] 


Cu(DC0 2 ) 2 -4D 2 


0.23(1) 




+ "V f\ ■ 

L ^> V 


iu 


z 




16.54(5) 


A 


103 


Rb 2 CrCl 4 


0.230(2) 


2 


in- 1 

1U 


j . 




-2 


52.3 


F 


[25] 


Gd 2 Cu0 4 


0.23 


7 


1 n-l 
1U 


3- 


10" 


-3 


6.4 


A 


104 


(C 6 H 5 CH 2 NH 3 ) 2 CrBr 4 1 


0.23 


7 


in 1 

iu — 


1- 


10- 


1 


52.0(1) 


F 


105 


KMnP0 4 -H 2 0t 


0.23(2) 




t > 9- 


10 


-2 




~ 15 


A 


93 


(CH 3 NH 3 ) 2 MnCl 4 


0.23(2) 


1 


10- 2 - 


1- 


10" 


-3 


44.75 


F 


106 


ErCl 3 


0.23(2) 


4 


10- 1 - 


1- 


10" 


2 


0.350(5) 


A 


107 


(d 6 -5CAP) 2 CuBr 4 


0.23(4) 


4 


lO- 2 - 


6- 


10" 


-3 


5.18(1) 


A 


108 


Li 2 VOSi0 4 § 


0.235(9) 


4 


10- 1 - 


2- 


10" 


2 


2.85 


A 


109 


Li 2 VOGe0 4 § 


0.236 












1.95 


A 


[TTO] 


Lao.o4Sr 2 .g6Mn 2 C>7 5 


0.24(2) 












145.0(5) 


A 


[m] 


Lao. 525 Sri. 475 Mn04 


0.24(3) 












110(1) 


A 


[TT2] 


RbFeF 4 


0.245(5) 


6 


10- 1 - 


1- 


10" 


-2 


133(2) 


A 


[90] 


MnPS 3 


0.25(1) 




t > 3- 


10 


-2 




78.6 


A 


[TT3l[TT4] 


Co 5 (OH) 8 (chdc) -4H 2 


0.25(3) 










-3 


60.5 


F 


[TT5] 


YBa 2 Cu 3 6 + I , 


0.26(1) 


5 


10- 2 - 


5- 


10" 


410 


A 


[TT6] 


Rb 2 MnCl 4 (B > 5.8 T) 


0.26(1) 


1 


10- 1 - 


2- 


10" 


-3 


54 


A 


[Hani] 


Rb 2 CrCl 3 Br 


0.26(1) 


9 


10- 1 - 


1- 


10" 


-2 


55 


F 


[TT71ITT81 


Rb 2 CrCl 2 Br 2 


0.26(1) 


9 


10- 1 - 


3- 


10" 


2 


57 


F 


[TT71ITT81 


KMnF 4 


0.26(1) 


3 


10- 1 - 


3- 


10" 


-2 


5.2(1) 


A 


[119] 


RbMnF 4 


0.26(1) 


3 


10- 1 - 


3- 


10" 


-2 


3.7(1) 


A 


119 



* Studied by Mossbauer spectroscopy. 
t Studied by bulk magnetometry. 

t Studied by proton Nuclear Magnetic Resonance (NMR). 
§ Studied by muon Spin Rotation (/iSR). 
■ Studied by ac susceptibility. 



Table A2: Summary of transition temperatures T c and magnetisation critical 
exponents /3 for epitaxial magnetic films grown on a range of substrates. The 
thickness d m ; n denotes the thickness at which these quantities were measured 
and t range indicates the range of reduced temperature over which j3 was 
measured. The magnetic anisotropy is indicated by the direction of the easy 
axis, and can either be perpendicular (_L) or parallel (||) to the film plane. 



System 


Structure 


d m in (ML) 


P 




t range 




T c (K) Anisotropy 


Method^ 


Reference 


Fe on 
Pd(100) 


bet, lxl 


2.0 


0.125(10) 




t < 3- 10~ 2 




613.1 


± 


ECS 


C2Q] 






1.2 


0.127(4) 


1 


io- 1 - 3-ur 


-3 


< 100 


± 


MOKE 




Ag(100) 


bec, lxl 


2.5-2.7| 


0.124(2) 


1 


lO" 1 -! ■ 10" 


-3 


324 


± 


MOKE 


[L221I29] 


W(110) 


bee, lxl 


0.8 


0.124(1) 


1 


10" 1 -4- 10" 


-3 


221.1(1) 


[no] 


SPLEED 


[123114] 






1.0 


0.134(3) 


1 


10 _1 -5 -10- 


-2 


224 


[110] 


SPLEED 


ED] 






1.7 


0.13(2) 








317 


[110] 


MOKE 


124 


Ag(lll) 


bec, lxl 


l.S 


0.139(6) 


1 


lO" 1 -! ■ 10" 


-3 


~ 450 


II 


MOKE 


125 






2.0 


0.130(3) 


1 


io -1 - i • io _ 


-3 


~ 450 


II 


MOKE 


[1251 


Cu(100) 


fct, 4x1 


~ 2.5J 


0.17(3) 


1 


10 -1 - 1-10" 


-2 


370 


-4 


MOKE 


[T261IT2T1IT281 


W(110)} 


bee, lxl 


0.82 


0.18(1) 


3 


10 _1 -1 ■ 10" 


-2 


282(3) 


[110] 


TOM, CEMS 


[1291 [1301 [131] 


Cu g4 Al 16 (100) 


fee, lxl 


4.0 


0.212(5) 


3 


10 -1 - 1 • 10" 


-2 


288(2) 


II 


LMDAD 


[1321 


W(100) 


bee, lxl 


1.6 


0.217(2) 


1 


10 -1 - 1- 10" 


-2 


207.8(1) 


[001] 


CEMS, SPLEED 


[321 IT33] 


Au(100) 


bee, lxl 


1.0 


0.22(1) 


1 


10 -1 - 1-10" 


-3 


300 


[001] 


SPLEED 


[134] 






2.0 


0.25(1) 


2 


10 _1 -1 ■ 10" 


-4 


290 


[001] 


ECS 


135 


W(100)J 


bee, lxl 


1.5 


0.22(2) 








282(1) 


[001] 


CEMS 


® 


V(001) 


bec 


3 


0.23(1) 


2 


10 -1 - 2-10 - 


2 


~ 190 




MOKE 


136 


Pdt. 




0.2-0.4}: 


0.23(1) 


2 


10^-2 -10" 


-2 


> 50 


II 


MOKE 





GaAs(lOO) 


bee, 2x6 


3.4 


0.26(2) 


1 


10" 1 -! ■ 10~ 


-3 


254.8(2) 


II 


MOKE 




Co on 
Cu(lll) 


fee, lxl 


1.0 


0.125 








433 


± 


TOM 


[T38] 






1.5 


0.15(8) 








460 


± 


MOKE 


[1391 


Cu(100) 


fee 


2.0 


0.24 








~ 240 


II 


MOKE 


[1401 [141] 



Continued on next page 



Table A2 — continued from previous page 



System 


Structure 


d min (ML) 


13 




t range 


Tc (K) 


Anisotropy 


Method} 


Reference 


Ni on 
W(110) 
Cu(lll) 
Cu(100) 


fee, 7x1 
fee, lxl 
fee, lxl 


2.0 
2.0-3.01 
4.1 


0.13(6) 
0.24(7) 
0.23(5) 


1 • 

3 
3- 


io -1 - i-io -3 

10- x -6- 10" 3 
10 -1 - 1-10 -2 


325 
380 
284 


II [001] 
II 

II 


FMR 
MOKE 
MOKE 


[142] 
[1431 
[Till 139 


V on 
Ag(100) 


bee, lxl 


5.0 


0.128(10) 


3- 


10 _1 -2-10- 4 


475.1 


II [001] 


ECS 


ma 


Gd on 
Y(0001) 


hep 


1.0 


0.23(5) 


1- 


10 _1 -8- 10" 3 


156 


II 


MOKE 


ma 


MnsGes on 
Ge(lll) 




1.0 


0.244 


2- 


10 _1 -4-10 -3 


296 


II 


SQUID 


[1471 


CoAl(lOO) 


bee, lxl 




0.22(2) 


2- 


10" i -7- 10~ 3 


~ 90 


II 


MOKE 


.2 



} Experimental properties were measured by the following techniques: Electron Capture Spectroscopy (ECS), Magneto Optical Kerr Effect (MOKE), Spin 

Polarised Low Energy Electron Diffraction (SPLEED), Torsion Oscillation Magnetometry (TOM), Convertion Electron Mossbauer Spectroscopy (CEMS), 

Linear Magnetic Dichroism in the Angular Distribution of photoelectron intensity (LMDAD) and Superconducting QUantum Interference Device (SQUID). 

} Exponent determined by averaging over values of a range of films of different thickness. 

} Reversible spin reorientation transition from || to _L with increasing T. 

} Coated with Ag. 

} Pd layers <5-doped with Fe. 
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Table A3. Chcmisorbed and physisorbed systems displaying two-dimensional 
phase transitions. 



System 




Model Ascribed 


Method 


Reference 


W(011)p(2 X 1)-H 


0.13(4) 


2d Ising* 


LEEDT 


GEO 


W(011)p(2 x 2)-H 


0.25(7) 


2dXY* 


LEED 


.36. 


p(l x 2) «(lx l)-Au(HO) 


0.13(2) 


2d Ising 


LEED 


A 


W(112)p(2 x l)-0 


0.13(1) 


2d Ising 


LEED 


m 


p(2 X 2)-0 disordering on Ru(OOl) 


0.13(2) 


3-state Potts 


LEED 


18 


p(2 X l)-0 disordering on Ru(OOl) 


0.085(15) 


4-state Potts 


LEED 


C3 


p(2 X 2)-0 disordering on Mo(110) 


0.19(2) 


XYhl 


LEED 


37 


p(2 X l)-0 disordering on W(110) 


0.19(5) 


XYhl 


LEED 


EU 


Ru(001)p(2 x 2)-S 


0.11(2) 


4-state Potts 


LEED 


33 


Ru(001)(\/3x V3)R30°-S 


0.14(3) 


3-state Potts 


LEED 


33 


(3 X 3)-Sn disordering on Ge(lll) 


0.11(1) 


3-state Potts 


HAS*, XRD§ 


34 


(3 X 1) reconstruction on Si(l 13) 


0.11(4) 


3-state Potts 


LEED 


m 


Xe melting on graphite 


0.23(4) 


2dXY 


XRD 


m 



* Model not ascribed by original authors, 
t LEED: Low Energy Electron Diffraction, 
t HAS: Helium diffraction. 
§ XRD: X-ray diffraction. 
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